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Abstract 



The capacitance of the circular parahel plate capacitor is calculated by expand- 
ing the solution to the Love integral equation into a Fourier cosine series. Previously, 
this kind of expansion has been carried out numerically, resulting in accuracy prob- 
lems at small plate separations. We show that this bottleneck can be alleviated, 
by calculating all expansion integrals analytically in terms of the Sine and Cosine 
integrals. Hence, we can, in the approximation of the kernel, use considerably larger 
matrices, resulting in improved numerical accuracy for the capacitance. In order to 
improve the accuracy at the smallest separations, we develop a heuristic extrapo- 
lation scheme that takes into account the convergence properties of the algorithm. 
Our results are compared with other numerical results from the literature and with 
the Kirchhoff result. Error estimates are presented, from which we conclude that 
our results is a substantial improvement compared with earlier numerical results. 

1 Introduction 

The exact capacitance of the circular parallel plate capacitor, with infinitely thin plates, 
remains an unsolved problem in potential theory, in the sense that to this date no explicit 
analytical solution has been reported. However, the problem can be formulated as a 
Fredholm integral equation of the second kind, known as Love's integral equation [1], 
which can be solved numerically. 

To our knowledge, the up to date most accurate studies of the capacitance at small plate 
separations are the ones by Wintle and Kurylowicz [2] and by Carlson and Illman [3]. 
Both of these studies have been used as benchmarks for solutions obtained by other 
methods; see e.g. [4-6]. Wintle and Kurylowicz [2] use an El-Gendi method [7] to rewrite 
the Love equation and apply numerical integration using the Clenshow- Curtis quadrature 
method [8] to obtain the capacitance. Carlson and Illman [3] solve the Love equation 
through an expansion of the kernel into a Fourier-cosine series. Later, they have also 
extended that method to solve the three-plate problem by means of coupled Love type 
equations [9]. It is known [1] that for small plate separations, a solution obtained via 
a series expansion of the kernel converges slowly, requiring a large number of expansion 
terms. To calculate the expansion coefficients of the kernel, Carlson and Illman [3] use 
numerical integration. Hence, their method is limited by a combination of the accuracy of 
the integration and the large number of terms needed. The accumulated errors effectively 
limit the expansion to about 100 terms, which is insufficient for convergence at very small 
separations. 

In this paper, we show that all of the integrals in the series expansion in [3] can be ex- 
pressed analytically in terms of the well-studied Sine- and Cosine integrals. In this way, 
we improve the numerical accuracy of the expansion coefficients up to the accuracy of the 
evaluations of the Sine- and Cosine integrals, which makes it possible to increase the num- 
ber of expansion functions considerably. Hence, in our method, the numerical accuracy in 
the capacitance is mainly limited by the truncation of the number of expansion functions. 
Thus, we can present improved results for the capacitance at small plate separations; 
results that surpass the results in [2,3], both in correctness and in the significant numbers 
of digits. Our results are also in excellent agreement with the result by Kirchhoff [10], 
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which becomes increasingly accurate when the plate separation tends to zero [11,12]. 

The paper is organized as follows: In Section [21 we review the expansion method, originally 
presented in [3]. In Section [3l we derive the analytical expressions for all the expansion 
integrals. Numerical results are presented in Section |4] and Section [5] contains some 
conclusions. 



2 Problem formulation and initial analysis 

The circular parallel plate capacitor is depicted in Figure [H 

The distance between the circular plates is denoted d and ! 

their common radius is denoted a. The model is idealized in I » 

the sense that the plates have zero thicknesses. Following d \ 
the notation used in many of the previous studies of this 

problem, we let k = d/a denote the normalized separation j 
between the plates. 

, . ^ , 11 1 1 ^- • FQi Figure 1: The circular par- 

ihe capacitance of the parallel plate capacitor is 3 n i i ^ / • i 

allel plate capacitor (side 

ri view slightly from above). 

C = A6oa / f{s)ds, (1) 

where the function /(s) is the solution to the modified Love integral equation 

f{s) - C K{s, t) fit) dt = l, 0<s<l, (2) 



with kernel 



K{s,t) = - 

TT 



(3) 



Note that in the original derivation by Love [1], the kernel and the function are defined in 
the range — 1 < s, t < 1, and the kernel has only one term, but since f{s) can be shown to 
be even one can instead use the formulations ([2]) and ([3]) (an elegant and short derivation 
of Love's integral equation can be found in [12]). 

To solve Equation ([2]) numerically, we follow the approach in [3] and expand the kernel 
and the unknown function into the Fourier-cosine expansion functions 



V^m(s) = V2 - Smocos(m7rs) , m = 0,l,..., (4) 

which in our study have been normalized to fulfil the orthogonality relation 

1 

(s) (s) ds = 5 w , (5) 







where bmm' denotes the Kronecker delta function. Note that similarly to f{s) all ipm- 
functions have vanishing first derivatives at s = 0. 
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Carrying out the expansions of f{s) and K{s,t), in terms of {ipm}, we obtain 

oo 

fis) = fmlpmis) , (6) 

m=0 

fm= f{s)iJrn{s)ds, (7) 







oo oo 



m=0 n=0 

K^n= [ [ K{s,t)^Us)Mt)<isdt, (9) 



which yield the following infinite linear system of equations for the coefficients {fn}°° 



n=Q- 



- Kmn) fn = ^mO, 171 = 0,1,.... (10) 



n=0 



From ([T]), (j6]) and the orthogonality of (jlj), the capacitance reduces to 

C = Asoafo, (11) 



where /o is simply the (0,0)-element in the inverse of the matrix with elements bmn~K„ 

In the numerical implementation the matrix is truncated into the size < m,n < N, 
where we call the truncation number. It is well-known [3,11] that at small separations, 
K, large values of are required to obtain convergence of /o in ( ITT]) . Also, for the accuracy 
of the result, it is crucial that the matrix elements Kmn, given by the integrals in IQ, 
have been calculated with a high accuracy. 

For small values of k the kernel K{s, t) has a pronounced crest ai s = t, making it difficult 
to evaluate the integrals in (Q numerically. In [3] this problem has been alleviated, by 
adding and subtracting a suitable term to the kernel, thereby removing the crest in one 
of the integrals, and making the inner integral in the other integral available for explicit 
evaluation. A related procedure has also been used earlier in [2]. However, a remaining 
problem is that for large values of m, n the expansion functions (jlj) yield rapidly oscillating 
integrands, which (when integrated numerically) result in slow convergence and poor 
accuracy. Hence, it would be beneficial if all integrals, encountered when expanding the 
kernel, could be expressed analytically in terms of established functions, readily available 
for numerical evaluation with a high accuracy. In the next section we will show that this 
is indeed the case. 



3 Expansion of the kernel into Sine and Cosine inte- 
grals 

In this section, we derive the analytical expressions for the expansion of the kernel K{s, t). 
First, we notice from and that the expansion coefficients have the property 
Knm = Kmn, resulting in a symmetric matrix. For the special cases when the indices 
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coincide and/or becomes zero, it is advantageous for the numerical evaluation to derive 
special simplified expressions for Kmn- In niost of the derivations, we utilize the following 
properties of the Sine and Cosine integrals [13]: 



Si(^*) = [Si(^)]*, 
Ci(z*) = [Ci(z)]*, 
Si(— 2;) = — Si(2;) , 

Ci(-^) = Ci(^) - jTT (0 < arg^ < tt) , 
where * denotes complex conjugation and j denotes the imaginary unit. 



(12) 
(13) 
(14) 
(15) 



3.1 Case m = n = 

By elementary integrals, we obtain 

•1 rl 



0,0 



K{s, t)dsdt^- [ I [ 
JO ^ Jo VJo 

1 



+ 



1 

TT JO 
1 

2^ 



arctan 



t + l 



K 



arctan 



t-1 



K 



dt 



4 arctan — — k m 1 H — - 



ds \ dt 



(16) 



3.2 Cases m = 0,n > 

Here the inner integral is the same as in the previous case, resulting in 

Kon = V2[ [ K{s,t)cos{mTt)dsdt 
Jo Jo 

— 1 M 

cos(n7rt) dt. 



V2 

TT Jo 



(t+l\ (t-l 

arctan — arctan 

\ K J \ K 



Defining the function 



/i(/3, K,a) — j cos{(3t) arctan + 



it follows that 



V2 



Kon = — [h{mT, 1) - /i(n7r, -1)] 

TT 



(17) 



(18) 



(19) 
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By a change of variable and integration by parts, we obtain 
ii(p, fi;,a)= / cos(pt) arctani I dt = < m = 

= K / cos(/?(km — a)) arctan('u) dw 

sin(/?) + 1 /•(^+")/'^ sin(/?(KM-a)) 



(3 



fl + a\ 1 ^^+"^/" sin a)^ 
arctan ~ tt / 5 

V / P ia/^c 1 + 



arctan 



(^) 



ilm|sin(/?(a (^Ci(/3 (a + Jk)) - Ci[(3{a + 1 

cos(/?(a+j/€))(^Si(/3(a+jK)) - Si(/5 (a + 1 + j/t))) | . (20) 



The evaluation of the last integral was carried out in the Maple software, and the result 
was simplified using the trigonometric-hyberbolic addition formulas and the properties 
([H and ([H. 



Insertion of fl20l) into f|T9l) and further simplifications yield 
V2 



mTTK] 



Kon = — 2-'"'^| cos(n7r (1 + jK))Si(n7r (2 + cos(n7r (1 — jK))Si(— j? 

— sin(r;,7r (1 +jK))Ci(n7r (2 sin(r?,7r (1 —jK))Ci(—jr27rK) k (21) 



3.3 Cases m n,m > 0,n > 

Here, we obtain 

Kmn = 2 / / -ft'(s, t) cos{mrt) cos(m7rs) dt ds 
Jo Jo 

cos(n7rt) cos(m7rs) dt ds. (22) 



^0 
1 /•! 



2 

JO JO 

Defining the function 



K K 
+ 



-r , X KCOS(3t) , . 

l2i/3,K,a)= ——^^dt, (23) 
Jo ft"' + (t + a) 

it follows that 

2 /"^ 

Kmn = - [hinn, k, -s) + hinn, k, s)] cos^mns) ds. (24) 
Jo 

Again, using Maple and some simplifications it follows that 

hiP, K, a) = Im |sin(/3 (a + Jk)) (si{p {a +]/«)) - Si(/5 (a + 1 + j/t))) 

+ cos(/3(a + j/s:))(^Ci(/5(a + j/€)) -Ci(/3(a + l+j/s:)))}. (25) 



Using dHD-dlSD, it follows that 

l2{P,K,-s) + l2iP,K,s) = -Imj sin(/5(s +jK))Si(/?(s + 1 

+ cos(/?(s +jK))Ci(/5(s + l+jn)) 
+ sin(/?(s-j«:))Si(/?(s-l-j«:)) 
+ cos(/? (s — jfi;))Ci(/5 (s — 1 — j/t 

Inserting f l26l) into the expression f l24l) for -ft'rrtn, we can write 

2 



vr 



(26) 



(27) 



where the function 

W, ^' ^) = ^ [^2(/5, -s) + /2(/3, s) ] cos(7s) ds = -^Imj 



sin((/5 + 7) s + j/t) + sin((/5 - 7) s + jfi:)) Si(/5 (s + 1 + jfi:)) 

+ (^cos((/3 + 7)s+jft) + cos((/5-7)s+jft))ci(/3(s + l+jft)) 

+ (sin((/3 + 7)s-j/€) +sin((/?-7)s-jfi:))Si(/5(s-l - jfi:)) 
((/5 + 7) s - j/t) + cos((/3 - 7) s - j/c))ci(/5 (s - 1 - j/c)) 



+ (^cos((/5 + 7) s - j/t) + cos((/3 - 7) s - Jk) jCi(/5(s - 1 - Jk)) ds 

= -^Ini{/4(/3 + 7,/3,jK, l+j/t) + /4(/5 + 7,/3, -j/?, -1 -j/t) 

+ hiP - 7,/3,jK, 1 + Jk) + h{P - -f,l3, -iK, -1 - Jk) }, 

and where 

/4(g, P, zi, Z2) = / [sin(gs + Pzi) Si(/?s + PZ2) + cos(gs + /?2;i) Ci(/3s + /5z2)] ds 
Jo 



(2^ 



sm{f3zi + q) Ci(/? (1 + Z2)) - cos{(3zi + q) Si(/5 (1 + Z2)) 
— sm.{[3zi) Ci{[3z2) + cos(/3zi) Si(/9z2) 
+ cos(/3zi - qz2) (Si((/? - g) (1 + ^2)) - Si((/3 - g) ^2] 
sin( ^ X ^- // 



in(/3^i - gz2) Ci((/? - g) (1 + ^2)) - Ci((/3 - q) Z2 



(29) 



To derive (1291) . we utilized the trigonometric addition formulas and the integration for- 
mulas 5.31 and 5.32 in [14]. Summarizing, we obtain 



Kmn = Im{/4((ra + m) 7r,n7r, Jk, 1 + Jk) + ^{{n + 



m] IT, nir, 



7T 



+ li^i^ri — m) 7!-,mT,}K, 1 + }k) + Ii[{n — m) TT,mT, — j/t, —1 (30) 



3.4 Cases m = n > 



In cases when m = n > the last two terms in (!30|l are not suitable for numerical evalu- 
ation, since in the final expression in (129|) the limit q ^ must be taken. Alternatively, 
we can start with setting g = in the integrand in (129|) : the result reduces to 

Knn = Im| l4{2nTr, nvr, Jk, 1 + Jk) + Ii{2mr^ mv, — Jk, —1 — }k) 

71 { 

+ sinQnTT/t) J ^Si(n7r (s + 1 + jfi;)) — Si(n7r (s — 1 — jfi;)) jds 
+ cosQnvr/t) y ^Ci(n7r (s + 1 + j/t)) — Ci(ri7r (s — 1 — j/t)) jdsj- 
= Im| l4{2mT, nvr, j/t, 1 + j/t) + /4(2n7r, nvr, — Jk, —1 — j/t) 



+ j sinh(n7r/t) (2 + Jk) Si(n7r (2 — j/tSi^nvrfi;) 

+ cosh(n7rK) ^ (2 + Jk) Ci(r27r (2 — jKCiQ^TT/t) — jvrj j-. (31) 



4 Numerical calculations 



The here presented improved expansion is most useful at small separations, wherefore we 
restrict our study to values k < 0.01; results obtained at larger separations are readily 
available in the literature, with some of the most accurate in [2,3]. 

The value of the truncation number is mainly dictated by the computational resources 
at hand. Here, we have used a maximum value of = 15000, at the smallest separations. 

All our results will be presented in terms of the normalized capacitance C = C/ {Aeoa), 
where Aeoa is the capacitance between two infinitely separated disks. Hence, it follows 
from ([n]) that C = /q. 

4.1 Extrapolation schemes 

To improve our numerical results, we have employed extrapolation. 



4.1.1 Power law model 

First, we considered a simple power law model for the capacitance: 

/o(Ar)=C-/5(Ar«:)-", (32) 

where fo{N) is the result obtained when using the truncation A^ in Equation (ITU]) . C is 
the extrapolated estimate in the limit N oo; a and /5 (assuming a > 0,/3 > 0) are 
coefficients that are determined together with C. 
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Nk = 0.2 


Nk = 1 




r> 




/L? ■ iU 










U.Ul 


80.235 


80.539 


1.380 


3.30 


O d O "1 o 

80.4312 


80.4363 


2.358 


5.15 


0.005 


ioo.yo i 


ioy./4U 


1 /I Q 1 
i.4oi 


o.Oo 


1 c;n 1 QQ/I 
ioy.iC)o4 


ioy. i4oo 


Z.o 1 1 


1 /I 

0. 14 


u.uuz 


394.778 


395.078 


1.472 


2.81 


394.9827 


394.9878 


2.384 


5.16 


0.001 


787.647 


787.946 


1.486 


2.73 


787.8533 


787.8585 


2.386 


5.16 


0.0005 


1573.217 


1573.516 


1.492 


2.71 


1573.4238 


1573.4290 


2.388 


5.16 


0.0002 


3929.640 


3929.849 


1.496 


2.69 


3929.8466 


3929.8518 


2.389 


5.16 


0.0001 


7856.804 


7857.102 


1.498 


2.69 


7857.0105 


7857.0156 


2.391 


5.15 



Table 1: The extrapolation parameters, a,j3 and C, obtained when fitting the power 
law extrapolation model (!32|) to fo{N) ,/o(A^/2) and fo{N/3). The estimate of the true 
capacitance is C = isoaC 



In our initial tests, the extrapolation parameters, a,P,C, were determined by fitting 
f l32|) to fo{N) ,/o(A^/2) and fo{N/3) (using rounded values for the fractions of when 
necessary). In Tabled] we present our results, for various separations n, for two different 
values of the product Nk. 

An important observation in Table [T] is that for constant values of Nk the values of the 
parameters a and f3 are approximately constant, regardless of the value of k. This is 
especially true for Nk = 1, but also for the smaller K-values when Nk = 0.2. The results 
for Nk = 0.2 (farther from convergence) and Nk = 1 (closer to convergence) together 
indicate that the extrapolation scheme over-estimates the capacitance. Another observa- 
tion from the data in Tabled] is that for each value of Nk the amount of extrapolation is 
approximately the same, regardless of the value of k. 

4.1.2 Heuristic model for improvement at low accuracies 

Considering decreasing separations k, the limitation of the truncation A^ will eventually 
force us to use smaller values of Nk, which will push us farther from convergence. However, 
the properties of the convergence, as demonstrated in Table dl offer an opportunity to 
improve poorly converged results, by using the following heuristic model: 

fo{N)=C + h{NK), (33) 

where the term —f3 {Nk)~'^ in (1321) has been replaced by a general function h{NK). Using 
Equation (1331) . we have developed a heuristic extrapolation method that is based on the 
assumption that when going from one K-value to the next smaller value the function h 
remains almost unchanged. To verify this, we have compared the /i-functions obtained at 
three different values of k. We used equation (|33]) . setting C = C, with C taken from the 
seventh column in Tabled! The results are given in Figure [21 as function of (A^k)~^ for the 
values between the ones in Table d] The /i-functions in Figure [2] essentially illustrate the 
convergence, with higher degree of convergence at lower {Nk)~^ values, where the /i-values 
at (A^k)~^ = 1 indicate the amount left to extrapolate. The three curves are essentially 
occupying the same distance, especially those at the smaller separations k = 0.001 and 
K = 0.0001, and for easier comparison we have included a magnification of the curves. 
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The nearly overlapping curves indicate that the convergence of the matrix inversion is 
determined mainly by the product Nk. 



-0.05 



h{NK) 




-0.15 



Figure 2: Comparison of the estimated /i-functions, obtained from Equation (133|) . for 
three different orders of the separation k. The /i-function essentially indicates how close 
to convergence the calculation of C is, flatness of the curve for high Nk values indicates 
that C has essentially the distance \h{Nk)\ left to its approximately true value. 

Now, having confirmed the assumption, the heuristic extrapolation method works as 
follows: 



1. Order the considered k- values as kq > ki > K2 > • • •, with N, fo{N) ,C and h 
indexed analogously. 

2. Starting with hq, for which we already can assume a fairly good convergence by using 
a moderate number of expansion functions (see Figure [2]), and using the power law 
f l32|) . with a truncation number Nq, we obtain the extrapolation Cq. This is taken as 
the starting value, Cq = Cq, for our improved algorithm, which from now proceeds 
repetitively: 

3. Increment the index (denoted i). Use equation fl33|) to find the /i-function from the 
previous step: 



Let ni_i be a truncation number fulfilling 

NiKi = ni_iKi_i. 
Assuming that hi = /ij-i, it follows from fl33|) and fl35l) that 



foA^i) ~ — fo,i-l(,^i-l) ~ ^i-l — /o,i-l 



Ni -C,_i, 



(34) 



(35) 



(36) 
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from which the extrapolated capacitance becomes 

C^ = hANi) + CVi - /o,^-l i^N^ ■ (37) 

Note that nj_i, determined from fl35l) . must fulfil ni_i < Ni_i, and if ni_i does not 
becomes an integer the last term in (IS7|) must be evaluated by interpolation between 
the adjacent integer values. 

4. Repeat step 3 until the final (smallest) K-value has been considered. 
4.2 Results for the capacitance 

At small separations, the first approximation to the capacitance is the geometric capaci- 
tance Cg = IT /{An). A much better approximation is the result by Kirchhoff [10]: 

^ £ + i ^^(^) + i t^^^^^'^) - ^] + • ^2^) 

This formula, which has been proved rigourously by Hutson [11], becomes increasingly 
accurate as k decreases. 

Our results are divided into two cases. In the first case, we used a constant value Nk = 3, 
and considered separations down to k = 0.0002; the results are given in Table [2l Here, 
all extrapolations, C, were obtained by fitting the power law (!32|) to the /o-values at 
N, N/2, N/3. We also present the relative excess (C — Cg)/Cg over the geometric capac- 
itance. Our results are compared with the Kirchhoff result ( l38l) and with the numerical 
results in [2, 3]. 



n 


N 


/o(iV) 


C 


Exc. Cg 


Eq.dMD 


Ref. [2] 


Ref. [3] 


0.01 


300 


80.43440 


80.43451 


2.41 % 


80.42044 


80.4342 


80.43 


0.005 


600 


159.14169 


159.14179 


1.31 % 


159.13354 


159.13 


159.1 


0.002 


1500 


394.98596 


394.98607 


5.82 %o 


394.98206 


394.87 


395 


0.001 


3000 


787.85661 


787.85672 


3.13 %o 


787.85443 


787.6 


787 


0.0005 


6000 


1573.42707 


1573.42718 


1.67 %o 


1573.42588 


1573 




0.0002 


15000 


3929.84994 


3929.85005 


7.28 %oo 


3929.84944 


3928.9 





Table 2: Numerically calculated capacitances, where fo{N) is the non-extrapolated value 
and C is the extrapolant. The estimated true capacitance is C = AeoaC. 



In the second case, we considered even smaller separations, down to k = 0.00001; the re- 
sults are presented in Table [31 Since limited computer memory enforced a maximum trun- 
cation number N = 15000, all extrapolations were obtained with our heuristic method. 
As the starting value for the repetitive extrapolation, we took the C-value obtained at the 
smallest separation considered in Table [2], i.e. kq = 0.0002, Co = 3929.85005; see Section 

Km 

At the largest separation, k, = 0.01, our result agrees better with the reference numerical 
results than with fl38|) . but for the smaller separations our results are more close to fl38l) . 
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J U V / 


c 


Exc. Cg 


Eg. (1381) 


Ref. [2] 


0.0001 


15000 


7857.01294 


7857.01378 


3.86 %oo 


7857.01355 


7855.9 


0.00005 


15000 


15711.16055 


15711.16855 


2.04 %oo 


15711.16847 




0.00002 


15000 


39273.25402 


39273.34241 


0.87 %oo 


39273.34244 




0.00001 


15000 


78543.05664 


78543.42381 


0.46 %oo 


78543.42390 





Table 3: Numerically calculated capacitances, where fo{N) is the non-extrapolated value 
and C is the extrapolant. 



In Figure [31 we have plotted the difference between our results and the Kirchhoff result, 
which in effect is our numerical approximation of the rest term in Equation (!38|) . For 
0.0002 < K < 0.01 the difference is positive and decreases with k with the approximate 
behavior oc k^'^, but when reaching the smallest k- values the difference turns negative and 
increases in magnitude. This cannot by itself be taken as an error, since (!38l) is not the 
exact solution, neither a lower bound, but a solution that becomes increasingly accurate 
as K decreases. 



c-a 



10 



10 



10 



10 



10 



K 



10 



10 



Figure 3: The magnitude of the difference between our extrapolated results and the 
Kirchhoff results from fl38|) . as function of the relative separation k. Red and blue colors 
denote positive and negative differences, respectively. The circles are the results from 
Table [2] and the dots are the results from Table O 



We should also mention the result by Ignatowsky [2,15], which differs from Kirchhoff 's 
result in that the constant term is replaced by (In 8 — 1/2) /4, thereby falling below the 
Kirchhoff result by the amount of approximately 0.33447. Polya and Szeg [16] have shown 
that Ignatowsky's result is a sharp lower bound for the capacitance. Thus, we see in Table 
|3] that at K = 0.00001 our non-extrapolated result is below this lower bound, which is 
an indication that convergence has not been reached due to an insufficient number of 
expansion functions. However, the extrapolated result is above the sharp lower bound. 
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4.3 On the accuracy of the results 

Our numerical simulations have shown that for a fixed k the unextrapolated capacitance 
fo{N) increases with A^, and from Tables [1] and [2] it appears that when using the power 
law extrapolation formula (l32l) the extrapolated value C decreases with A^. If we in Table 
[H use the C-values at Nk = 1 as references, the extrapolations obtained at Nk = 0.2 
overestimate the references with about one third of the total amount of extrapolation. 
Similarly, if we in Table [2] use the C-values at Nk = 3 as references, the extrapolations 
obtained at Nk = 1 (in Table [1]) overestimate the references with about one third of the 
total amount of extrapolation. Applying this rule to the values in Table |2], we conclude 
that the true values are approximately 4 ■ lO'^ below the extrapolated value. 

Using the subsequent repeated extrapolation, the error obtained at k = 0.0002 is propa- 
gated to the lower K-values. If the extrapolation scheme was ideal, i.e. if the /i-function in 
( I33l) was the same regardless of k, no further contributions to the error would occur. To 
get a very rough estimate of the cumulated errors, we can study the similarities between 
the curves in Figure [2j Given the excellent agreement between the curves at the smaller 
K- values, a pessimistic estimate is that the cumulated error is within 1/10 of the amount 
of extrapolation. Applying this estimate to the values given in Table [3] we obtain the 
estimate of the maximum error, denoted AC, acquired in each step of the interpolation 
scheme. The results, rounded upward to one significant digit, are presented in Table HI 
Since we typically have to let Nk decrease with k, the most significant contribution to 
the accumulated error at a certain k- value is acquired in the last step of the extrapolation 
scheme. 





0.0002 


0.0001 


0.00005 


0.00002 


0.00001 


AC 


4 ■ 10-5 


8 ■ 10-5 


8 ■ 10-^ 


9 • 10-3 


4- 10-2 



Table 4: Estimated order of the cumulated error during each step of the heuristic ex- 
trapolation scheme. 



5 Conclusions 

By expanding the kernel in the Love equation into a Fourier cosine series, with the co- 
efficients expressed analytically in terms of Sine and Cosine integrals, we have increased 
the accuracy of the expansion, making it possible to use considerably larger truncation 
numbers N than in previous studies. In this way, we have improved the numerical values 
for the capacitance at small plate separations k. 

The present method enable us to consider smaller plate separation distances than has been 
considered before, while maintaining a high accuracy. Numerical tests indicated that the 
degree of convergence is determined by the product Nk. Hence, for this method of calcu- 
lating the capacitance it becomes practically impossible to accommodate for a decreasing 
plate separation by a corresponding increase of the number of expansion functions; cf. 
Figure [2l To compensate for this problem, we have developed a heuristic extrapolation 
scheme that uses the information obtained at intermediate separations, where the conver- 
gence is good, to improve the convergence at small separations. 
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At larger separations, k, than those considered in our study, convergence is obtained for 
smaller numbers, A^, of expansion functions. In such cases, one does not need to use 
the analytical results in Section [3] for the integrals, since they can instead be calculated 
accurately by means of numerical methods [2,3]. In fact, for large k- values the explicit 
expressions in Section [3] are unsuitable for numerical evaluation. This can be seen by 
observing that from and (jll) it follows that for any k > the matrix elements fulfil 
the relation 

\Kmn\< 2Kofl < 2 (39) 

and that when we tested the algorithm for separations 1 < k < 10 we encountered spurious 
large matrix elements violating the condition (139|) . The reason is that when Nk ^ 1, 
the individual terms in the expressions given in Section [3] exhibit an exponential growth 
oc e^*^ and the condition ( |39l) is met by taking the differences between such very large 
terms, but numerically this kind of evaluation leads to cancelation effects. 
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